void ENER(En, Vir, R2)
{
	double R2, r2i, r6i, En, Vir;
	if (R2 < RC2)
	{
		r2i = SIG2/R2;
		r6i = r2i*r2i*r2i;
		if (SHIFT)
			En = EPS4*(r6i*r6i-r6i) - ECUT;
		else
			En = EPS4*(r6i*r6i-r6i);
		Vir = ESP48(r6i*r6i-r6i);
	}
	else
	{
		En = 0;
		Vir = 0;
	}
}
